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We show that condensation in a capped capillary slit is a continuous interfacial critical phe- 
nomenon, related intimately to several other surface phase transitions. In three dimensions (3d), 
the adsorption and desorption branches correspond to the unbinding of the meniscus from the cap 
and opening, respectively and are equivalent to 2d-like complete-wetting transitions. For disper- 
sion forces, the singularities on the two branches are distinct, owing to the different interplay of 
geometry and intermolecular forces. In 2d we establish precise connection, or covariance, with 2d 
critical- wetting and wedge- filling transitions, i.e. we establish that certain interfacial properties in 
very different geometries are identical. Our predictions of universal scaling and covariance in finite 
capillaries are supported by extensive Ising model simulation studies in 2d and 3d. 



Capillary condensation (CC) is central to our under- 
standing of confined fluids and has received much atten- 
tion over the last few decades As is well known, 
vapour confined in a slit of width L condenses at a pres- 
sure p co (L), less than the value p sa t at bulk saturation, 
given by the macroscopic Kelvin equation 0, Q j 



Psat -Pco(L) = 



2<j cos ( 



(1) 



where a is the liquid-vapour surface tension and 9 the 
contact angle. Studies of CC, based on Landau or mod- 
ern density functional methods, usually consider confin- 
ing walls of infinite area and assume translational invari- 
ance parallel to these [3,0]. In this case, CC is certainly 
a first-order transition and mean-field treatments yield 
adsorption isotherms with a van der Waals loop. Far less 
attention has been given to CC in slits of finite depth D 
which are capped at one end and open into a reservoir 
(see Fig. 1). This scenario is certainly experimentally 
accessible and is similar to recent analysis of adsorption 
on grooved and pitted substrates 0, @ . Three numerical 
studies @, H, Q, restricted to complete wetting (9 = 0), 
have reported the same basic finding: In a finite capil- 
lary, CC is a sharp but continuous transition (Fig. 1 ) 
and adsorption isotherms exhibit no van der Waals loops 
for any finite D (see also [IS El)- However, the physical 
reason behind this striking change and the quantitative 
aspects of the transition have not been elucidated. 

Here we show that CC in such capped capillaries is 
a continuous interfacial critical phenomenon exhibiting 
observable critical singularities which are intimately re- 
lated to several other surface phase transitions. Our two 
main findings are: (A) In 3d, adsorption and desorp- 
tion in a deep capillary correspond to the continuous un- 
binding of the meniscus from the bottom and top, re- 
spectively, and map onto two-dimensional complete wet- 
ting [H, HI] with relevant scaling field Ap = p co (L) — p. 
The divergence of the average meniscus height on these 
branches is described by the respective critical singular- 
ities (£) ~ Ap~P A ,D - (I) - \ Ap\-P D which are, in gen- 



eral distinct. In particular, for dispersion forces, (3a = j 
while Pd = | implying that the adsorption (isotherm) 
is steeper than the desorption. For a finite-depth capil- 
lary, the meniscus has large-scale fluctuations at p = p co 
equivalent to those of an interface in a 2d infinite capil- 
lary with opposing walls for which there are long-standing 
predictions [14], HH HE Qj} • (B) In a 2d capped slit (or a 
3d capillary pore), one can proceed further and relate CC 
to critical wetting transitions occurring at 2d planar sub- 
strates [HI, [l8[ . This precise connection is an extension of 
the geometrical covariance known for the filling of wedges 
HiJ H3] , cones [2l| and apexes [22[ and implies that some 
(universal) interfacial properties in very different geome- 
tries are identical. These predictions, based primarily on 
analysis of interfacial Hamiltonians, are supported fully 
by our Ising model simulation studies in 2d and 3d. 
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FIG. 1: Schematic cross-section of a capped capillary of depth 
D and width L illustrating the local interfacial height ^(x). 
The slit is infinitely long in the y direction. A typical adsorp- 
tion isotherm is also sketched showing the rapid but contin- 
uous rise in the mid-point height near CC (p « p co (L)) and 
the complete wetting as saturation is approached. 

In an infinite capillary, the two phases, referred to 
as capillary-liquid and capillary- vapour, coexist when 
P = Pco(L). Any meniscus separating these phases is 
delocalised and has large-scale fluctuations, analogous 
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to those of a planar interface between coexisting bulk 
phases, but in one lower dimension. Now in a capped 
system, geometry necessitates the formation of a menis- 
cus whose location is determined, in part, by the capil- 
lary thermodynamics. For p < p co (L) the meniscus must 
be located near the bottom because the capillary-liquid 
phase is metastable, while for p > p co (L), the capillary- 
vapour is metastable and the meniscus must reside near 
the capillary opening. The fundamental question is, what 
happens to the meniscus as p — > p co l 

We first follow [H and introduce an interfacial model 
based on the local interfacial height £(x,y), where y is 
the variable along the capillary. The free energy is 

H[£] = aS[£] + S P V[£] + W[£] (2) 

Here S is the fluid interfacial area, V the volume of liquid, 
Sp = p sat —p, and W is the binding potential accounting 
for the (dispersion) intermolecular forces. For confining 
walls which are completely wet (9 = 0), one can use 

W[£] = ^ JJdrdv' (3) 

where A > is the Hamaker constant, an energy, and r 
and r' denote points in the substrate and vapour, re- 
spectively. The numerical prefactor is chosen so that 
([3]) recovers the usual binding potential per unit area 
W{£) = A£~ 2 for planar walls pjj. The model © 
describes the whole adsorption isotherm including the 
CC and the complete wetting of the entire substrate as 
V ~~ > Psat- The latter is not of interest so we seek an 
effective Hamiltonian which will allow us to study the 
CC more easily. Accordingly we integrate out degrees 
of freedom keeping only the long-wavelength fluctuations 
in the meniscus height along the capillary. Close to CC, 
the meniscus must be of semi-circular cross-section at lo- 
cal height £(y) = £(0, y) < D. The fluctuations of the 
meniscus are then described by a capillary Hamiltonian, 
obtained via a constrained minimisation of ©. After 
some algebra we find [23[ 

H cap [l(v)] = Jdy |^ ^ + W cap (£;D)^ (4) 

where the capillary binding potential is 

W cap (£;D) = A P L£ + + + ... (5) 

The model is only valid near CC and does not describe 
the complete wetting as p — > p sa t- Nevertheless it de- 
scribes accurately the interplay of the geometry, forces 
and fluctuations for p « p co , providing a more transpar- 
ent view of the interfacial behaviour. 

Eqn (HJ) is one of our central results and highlights 
some intriguing properties of CC: (I) The first term, con- 
jugate to £(y), is proportional to Ap = p co (L) — p, mea- 
suring the deviation from capillary coexistence. This es- 
tablishes the connection with complete wetting phenom- 
ena, in one lower dimension, with the meniscus acting as 



the unbinding interface. (II) Intermolecular forces re- 
pel the meniscus from the capillary cap and opening. 
At CC, the capillary cap wishes to be "wet" by the 
capillary- liquid while the capillary opening is "wet" by 
the capillary- vapour. The finite-size (FS) effects mimic 
therefore those of an interface in a capillary with opposing 
walls that are wet by different fluid phases 0, US 0, HtJ . 
This connection is deeper in 2d as we shall see later. 
(Ill) While the repulsion from the top is similar to the 
planar wall [251 ]. that from the bottom is shorter-ranged 
due to a geometry-induced cancellation of intermolecu- 
lar forces. This is the reason behind the asymmetry in 
the adsorption isotherm seen in numerical studies [1] and 
is a general feature of CC in systems with long-ranged 
(dispersion) forces. For short-ranged forces the (direct) 
repulsions from the top and bottom are similar to expo- 
nential decays. 

To continue, consider the critical behaviour in the 
semi- infinite limit D — > oo. The adsorption and des- 
orption branches of the 3d capillary isotherm become 
analogous to 2d-like complete wetting phase transitions 
describing the unbinding of the meniscus from the cap 
and open end, respectively. The associated critical sin- 
gularities are well understood [HI and reflect the long- 
ranged forces presented in W cap (£) and fluctuation-effects 
associated with the wandering of the meniscus controlled 
by the surface tension term in ((4|). The dependence on 
the slit width is significant since increasing L effectively 
supresses the role of fluctuations. For example, the en- 
tropic repulsion from the bottom effectively adds a term 
oc (kBT) 2 /aL£ 2 to ([5]). While, in principle, the asymp- 
totic divergence of £ is ultimately determined by this en- 
tropic repulsion, in practice it is irrelevant for slits more 
than a few angstroms in width because the amplitude is 
negligible compared to that of the long-range forces. The 
latter dominate for all practical purposes and we antici- 
pate mean- field-like behaviour with (£ } ~ (AL/ Ap) 1 / 4 , 
as quoted earlier. On the desorption branch, fluctua- 
tion effects are similarly negligible, even though they are 
marginal, and we predict D - (£) ~ (A/\ Ap]) 1 / 3 . For 
systems with short-ranged forces, on the other hand, the 
influence of fluctuations can no longer be neglected and 
both the adsorption and desorption branches show the 
same critical behaviour. Thus, for example on the ad- 
sorption branch (£} ~ (L 2 Ap) -1 / 3 . All these results are 
supported fully by transfer-matrix analysis of the Hamil- 
tonian (g]). 

Condensation (or rather pseudo-condensation) in a 2d 
capped capillary is particularly interesting. Here we fo- 
cus solely on systems with short-ranged forces since long- 
ranged forces are not pertinent for all physical scenar- 
ios. However, we broaden our discussion to the case of 
non-zero contact angle 9. Firstly, as above, the adsorp- 
tion and desorption in a semi-infinite capillary are related 
to the Id limit of complete wetting. The probability of 
finding the meniscus at a given height is simply propor- 
tional to the Boltzmann weight of W cap , implying that 
on the adsorption branch (I) ~ (LAp) -1 and similarly 
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for desorption. Thus the critical exponent 0a. = 1 can 
be identified as the d — ► 1 + limit of the general result 
Pco(d) — (3 — d)/(d+l) appropriate to short-ranged com- 
plete wetting in dimension d < 3 

To see the deep connection with 2d critical wetting we 
introduce what at first appears to be an artificial geom- 
etry, and adopt a magnetic notation also useful for com- 
parison with our Ising model simulation results. Consider 
a planar Ising-like system, with spontaneous magnetisa- 
tion mo and a boundary (wall) in the shape of a trough. 
That is the height of the wall above some reference line 
is %(x) = for |x| < L/2 and = (|a;| - L/2) tana 

otherwise. The spins on the boundary are subject to a 
surface field hi > 0. The spins away from the wall are 
subject to a position dependent external field h(x) = h 
for \x\ < L/2 and h(x) — Q~ otherwise. In general 
these boundary conditions induce a fluctuating interface 
at height £(x), which defines the "wetting" layer of up 
spins adsorbed near the wall. 

The trough has limiting geometries, each having dis- 
tinct transitions. For a = we have a capped capillary 
with spins subject to bulk field h exhibiting continuous 
CC when h = h co (L) « —a cos 8/(moL). In the limit 
L — > we recover a wedge, which at bulk two-phase co- 
existence exhibits a 2d filling transition [TiJ [20, [H, Hi! 
when 8 (T fin) = a, whereby the thickness of the adsorbed 
layer at the bottom diverges: l w ~ (Tfm—T)~P w . Finally 
if a, L — > we recover a planar geometry with a critical 
transition at temperature T wet , at which the mean layer 
thickness diverges: ~ (T wet — T)~P B , and 9^0 [12 ] . 

The trough exhibits a generalised filling transition 
which can be studied using the effective Hamiltonian 



H TR [l] = J dx{a 



1 ( M 

\dx 



W tr (i,x;h,a)] (6) 



The binding potential for the trough, Wt r , has a contact 
(attractive) interaction at the wall and a term 2m§h£(x) 
for | a; | < which models the bulk- like field acting on 
spins in the strip vertically above the bottom. The model 
is amenable to transfer matrix analysis similar to that 
used for filling in acute wedges [20j . In particular, one 
can establish that, for all a < n/2, close to the transition 
the probability distribution function (PDF) for the mid- 
point (x — 0) interfacial height has the simple scaling 
form Pt r {£) <x e^'^' where the mean- height is given by 



2(crsin(# — a) — mohL) 



(7) 



This result identifies correctly the phase boundaries for 
CC, filling and wetting in the limits discussed above. 
When L = 0, we recover the covariance between the 
PDFs at 2d wedge filling and planar critical wetting writ- 
ten, in an obvious notation, P w (£;9,a) = P„(£;6 — a) 
[l9| . The covariance of the PDF therefore extends to 2d 
CC and establishes a dimensional reduction between 2d 
critical- and Id complete-wetting for thermal systems. 
This has a number of consequences. First, the 2d expo- 
nents satisfy Pa = Pw — fl s — 1 and are identified with 



the Id complete wetting result /3* D (1 + ) noted earlier. Sec- 
ond, the invariance of the PDF necessarily implies that 
the thermal interfacial wandering exponent C(d) must 
satisfy C(2) = 1/2 and ((1+) = 1. These last results 
are not new but they point to a very deep connection 
between substrate geometry and interfacial fluctuations 




FIG. 2: Representative Monte Carlo simulation results 
(points) for the mid-point magnetization profile in 2d and 
3d capped capillaries compared with theoretical predictions 
©, © (curves). The temperatures, T = 0.741T c 6u? * (2d) and 
T = 0.8301^""% respectively, were chosen to avoid the bulk 
critical region and also the roughening transition (in 3d). 



Our discussion has focused so far on the critical sin- 
gularities in a semi-infinite capillary. Fluctuation effects 
and covariances are also apparent in a capillary of finite- 
depth, particularly if the pressure (or bulk field for Ising 
systems) is tuned to the condensation value, p = p co (L), 
of the infinite capillary. Here we focus on systems with 
short-ranged forces where the geometrical covariances are 
most striking and we can compare with simulation stud- 
ies. We consider only the case of strongly adsorbing walls 
(6 = 0), although the 2d case is insensitive to 9 and show 
that the fluctuations of the meniscus in both 2d and 3d 
capillaries are related, via geometrical covariance, to uni- 
versal FS effects in an entirely different 2d geometry. 

Consider an infinitely long, 2d Ising magnet with 
boundaries at z = 0, D along which the s pins are sub- 
ject to fields hi and —hi respectively [Til. |24|. In zero 
bulk field, the lower wall is wet by up spins for T > T we t 
while the upper is wet by down spins. The FS effects 
are known to fall into two universality classes. Exactly 
at the critical wetting temperature T = T wet the profile 
contains a scaling contribution [14| 



m(z) 
mo 



1 
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while for T 



bulk 



>T>T U 



2nz 



m(z) 2z 1 

— — = 1 h-sm. 

mo D 7T V D 



(8) 



(9) 
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These results are valid in the scaling limits z, D — > oo 
with z / D arbitrary and are verified by exact Ising calcu- 
lations [H, [13] ■ Both are indicative of large-scale inter- 
facial fluctuations whereby the roughness scales with the 
separation D. 

Now consider the magnetisation profile m(z) mea- 
sured along the central axis x — of a finite-depth 
capped capillary at bulk field tuned to condensation, i.e. 
h = h co (L). If the capillary is 2d, covariance with 2d 
critical wetting implies we should expect (jHJ). In a 3d 
capillary on the other hand, transfer matrix analysis of 
(HI, with short-ranged forces, leads us to the prediction 
©. 

In order to test these predictions Monte Carlo simu- 
lations were performed using a Metropolis algorithm for 
system sizes L = 13, D = 200 (2d) and L = 11, D = 88, 
M = 9000 (3d). The length M in the direction along 
the 3d capillary is sufficiently large to preclude FS. More 
quantitatively, transfer matrix analysis shows these are 

negligible if exp ^— f^gjf J << 1 where u) is the dimen- 

sionless wetting parameter [12[ and k the inverse bulk 
correlation length. We mimic the capillary cap geometry 
by fixing all spins on all surfaces to +1 except along the 
top line/plane z — D where they are fixed to —1. Peri- 
odic boundary conditions apply (in 3d) along the capil- 
lary. First we determined the CC line h co (L) for the 2d 
and 3d open capillaries, where periodic boundary condi- 
tions apply at the top (z = L) and bottom (z = 0), using 
standard multicanonical and histogram reweighting tech- 
niques. The measured form of the axial magnetisation 
profiles at h co are shown in Fig. [2j The comparison with 
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